KEK-TH-733 
KOBE-TH-00-10 



Pseudo-Dirac Scenario 
for 

Neutrino Oscillations 



Makoto Kobayashi^ and C.S. Lim^ 



( a ) Institute for Particle and Nuclear Studies, KEK, Oho 1-1, Tsukuba 305-0801, Japan 
^ Department of Physics, Kobe University, Nada, Kobe 657-8501, Japan 



Abstract 

We argue how pseudo-Dirac scenario for neutrinos leads to rich neutrino oscillation phe- 
nomena, including oscillation inside each generation. The pseudo-Dirac scenario is general- 
ized by incorporating generation mixings and formulae for the various neutrino oscillations 
are derived. As the application we compare the formulae with the corresponding data. We 
find that observed pattern of mixings, such as almost maximal mixing in the atmospheric 
neutrino oscillation, is naturally explained in the generalized Pseudo-Dirac scenario with 
small generation mixings. We, however, also point out that there remain some problems to 
be settled for this scenario to be viable. The possible theoretical framework to realize the 
pseudo-Dirac scenario is also briefly commented on. 



1. Introduction 

The presence of neutrino oscillation, strongly suggested by the recent Superkamiokande 
results on atmospheric neutrinos [|IJ, is almost unique clue to physics beyond the well- 
established standard model. More precisely, recent data on neutrino oscillations seem to 
have put forward the following challenging theoretical problems, which may lead to the 
physics beyond the standard model; 

(a) The data from (Super-)Kamiokande on atmospheric neutrinos necessitate a large or al- 
most maximal mixing angle Solar neutrino deficit may also be explained by the "large 
angle solution" |2| . How can such large or maximal mixing be naturally derived theoretically 

? 

(b) If we further accept LSND result ||, in addition to the solar and atmospheric neutrino 
data, the scheme with only 3 light neutrino states clearly gets into trouble. What kind of 
theoretical framework or model is needed to accommodate all of these neutrino oscillations 
? It has been argued that we should introduce at least one "sterile state" . 

(c) There seems to be a large disparity among the magnitudes of mixing angles implied by 
these experiments; The data on atmospheric neutrino and possibly that on solar neutrino 
indicate the necessity of large or almost maximal mixing angles |]J, ||, while the neutrino 
oscillation at LSND experiment is well described by a small mixing angle ||. How can such 
disparity be naturally explained theoretically ? 

One clear thing is that these problems altogether suggest that flavour mixing or mass 
matrices in leptonic sector are quite different from those in quark sector. In particular 
the presence of large mixings and the necessity of extending the scheme with only 3 light 
neutrinos are specific new features in the leptonic sector, not shared by the quark sector, 
and may lead to a drastic modification of the Standard Model. 

The main purpose of the present work is to generalize the pseudo-Dirac scenario for neu- 
trino masses, which we argue to imply rich neutrino oscillation phenomena, by incorporating 
generation mixings and to derive general formulae for neutrino oscillations in the generalized 
pseudo-Dirac scenario. We will show that the above problems, (a), (b) and (c), are naturally 
(without any fine-tuning) solved simultaneously in the generalized scheme without enlarging 
the number of generations. It is worth noticing that small generation mixings are shown to 
be exactly what we need to solve the problems. We, however, will also point out some serious 
problems encountered by the atmospheric and solar neutrino oscillations into sterile states 
[|J, @ and try to discuss some possible ways to cure these problems. We also investigate 
briefly what kind of model or theoretical framework is possible for the pseudo-Dirac scenario 
to be realized, for the purpose of searching for some direction to the physics beyond the stan- 
dard model. The key ingredient for our scenario is "pseudo-Dirac" property of neutrinos, 
whose precise meaning is now discussed in some detail. 

It will be natural to suspect that the specific features of flavour mixing or mass matrices 
in leptonic sector stated above should be related with the peculiarity of the leptonic masses, 
i.e. the fact that only neutrinos may have Majorana masses. Once neutrinos are allowed to 
have Majorana masses, we may think of three typical cases for neutrino masses, which we 
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will discuss successively below. In the base of weak-eigenstates, ip W L, where active states are 
put in the up-stairs and "sterile" states are put in the down-stairs, 

ip wL = ( - aL ) [a = e, /I, t; V aL = (v a R) C ), (1) 
the neutrino mass term is generally written as 

£mass = l -^ wL CMi) wL , (2) 
where C is the charge-conjugation matrix and the 6x6 mass matrix M takes a form of 

(3) 



\M D M* R 

with the 3x3 matrices Md, Ml and Mr being those for Dirac masses, left- and right- 
handed Majorana masses, respectively. Depending on the extent of lepton number violation, 
or relative magnitudes of Majorana masses to those of Dirac masses, we can think of three 
typical cases. 

(1) pure Dirac 

Imposing lepton number conservation, or ignoring all Majorana masses (Ml = Mr = 0), 
we get pure Dirac neutrinos. We all know that there are only 3 mass-eigenvalues for 3 
generations, although the mass matrix M should have 6 eigenvalues, in general. What's 
really happening is that there are 3 degenerate pairs of mass-eigenstates. The situation may 
be easily seen for simplified 1 generation case 

M=(° 7). (4) 

It is easy to see that M has eigenvalues mr>, —n^D an d the angle, 9, in the orthogonal matrix 
to diagonalize M is just 7r/4. If processes we are interested in have no chirality flip, as in the 
case of neutrino oscillations (in the absence of magnetic field), the sign of mass is irrelevant 
and we have degenerate mass-squared. It should also be noted that in this case a maximal 
mixing, 9 = 7r/4, between an active state and a sterile state, v L and Vl, has been realized. 
Unfortunately, this maximal mixing does not lead to any neutrino oscillation, just because 
the mass-squared are degenerated. 

(2) "Pseudo-Dirac" 

What happens if we allow small lepton number violation, i.e. if we switch on very small 
Majorana masses, Ml-, Mr <C Mo (with the magnitudes of matrices being compared by 
taking typical orders of magnitudes of the matrix elements). Neutrinos are still almost 
Dirac particles and are called "pseudo-Dirac" neutrinos [0. The small Majorana masses, 
however, slightly lift the degeneracy of mass-eigenvalues, and we get almost degenerate pairs 
of eigenstates with tiny mass differences. As far as the Majorana masses are small, the 
mixing angles should remain almost maximal, 9 ~ 7r/4. To understand the situation, we 
consider the 1 generation case again. The mass matrix now reads as 

M=( mL mD ), (5) 
\m D m R J 
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with ttil, m R <C mo- We now get |tan2#| = | rr ^ rr ^ n | 3> 1, leading to almost maximal mixing 
9 ~ ~. We have two mass eigenstates, which are almost symmetric and anti-symmetric 
combinations of active and sterile states, i.e., v$ = sin9 + cos9 T>l ~ 7f + an d 
va = (—i)(cos9 ul — sin# V L ) ~ ^ i ( u L — Vi). Their masses are almost degenerate but are 
slightly different; iris ~ m A ~ ffc, Am = 777.5 — in a ~ % + jtir <C mp. Now the tiny 
mass difference and the almost maximal mixing will lead to a neutrino oscillation between 
an active state and a sterile state, even if we have only f generation ! Actually even for the 1 
generation case two kinds of neutrino oscillations, without and with chirality flip, are possible 
(though the latter oscillation necessitates the presence of magnetic field and we ignore the 
magnetic field unless otherwise stated@in this paper), i.e. (i) v L — > v L and (ii) v L — ► v R ||. 
If matter effects are included, the oscillations of (i) and (ii) become resonant oscillations, 
which are quite similar to those in MSW J7| and Resonant Spin Flavour Precession (RSFP) 
|| scenarios: The Hamiltonian in the base of for the former resonant oscillation is 

given by @ 

H-( a ^r sin2e ) (6 ) 

H -[^sin29 ^cos2tfJ' (6) 

where the matter effect a is given for e.g. v e as a e = ^|(2iV e — N n ), with N e and N n denoting 
the number densities of electron and neutron, respectively. 

(3) See-saw 

The last possibility is famous see-saw scenario in which SU(2) invariant Majorana masses, 
M R , are supposed to be much larger than the Dirac masses: M R 3> M D , M L ~ 0. The sterile 
states Vl approximately become mass-eigenstates and are decoupled from low energy pro- 
cesses such as neutrino oscillation. Thus only lighter 3 mass-eigenstates (~ z/ a i) participate 
in neutrino oscillation phenomena. In the see-saw scenario, therefore, the mixings relevant 
for the neutrino oscillations are generation mixings and there seems to be no immediate rea- 
son to expect large mixing angles. It is also worth noting that as far as chirality preserving 
oscillations are concerned, there is no observable distinction between the cases of (1) and 
(3); In both cases only 3 light neutrino states participate in the oscillations. 

From the above discussion we learn that only in the pseudo-Dirac scenario 6 neutrino 
states fully participate in low energy processes, and rich neutrino oscillation phenomena, 
both inter-generational and active <-> sterile, are expected. In fact a trial to explain existing 
data on neutrino oscillations based on the pseudo-Dirac scenario was made some time ago 
||, fllQl . (For the recent revived interest in this scenario, refer to [jTT|| , [Q. See also Ref . [|13|1 
for the recent discussions in the framework of four neutrinos with only one sterile state.) 
In our previous attempt |§ generation mixings were switched off, for brevity, and three 
generations shared their roles to account for neutrino oscillations; Solar neutrino oscillation 
was mainly due to u eL — > V eL in the first generation (in ref. the effect of magnetic field was 
also taken into account). The atmospheric neutrino oscillation could be naturally explained 
by VfiL ~ > ^Vl with almost maximal mixing angle. At that time there was a datum to suggest 



the existence of 17KeV neutrino |14| , which enforced us to rely on the pseudo-Dirac scenario, 
as otherwise three generation scheme could not explain all of these data simultaneously. The 
17 KeV neutrino has been ruled out, and instead there has appeared the LSND data B, 
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which again necessitates three independent mass differences and to modify the ordinary 3 
generation scheme with either see-saw or pure Dirac neutrinos. 

In the present paper we generalize our previous argument || including generation mix- 
ings. We will find that small generation mixings, just as in the CKM matrix in the quark 
sector, is exactly what we need: the nice features of pseudo-Dirac scenario, such as maximal 
mixings in atmospheric and solar neutrino oscillations, are known to remain basically intact, 
while LSND data is naturally explained by the small generation mixings. In such a sense, 
we may say that 

"The recent data on neutrino oscillations may be natural consequences of the property that 
lepton number is only slightly violated and generation mixings are small. " 



2. Mass eigenstates and mixings 



In the present work we allow arbitrary generation mixings or arbitrary off-diagonal mass 
matrices, except for the assumption of pseudo-Dirac property M L ,M R <C M D . Thus we 
may naively expect that diagonalization of 6 x 6 mass matrix is quite complicated and 
various formulae for the probabilities of neutrino oscillations are expressed by use of 6 mass- 
eigenvalues and 15 mixing angles (together with possible multiple CP violating phases). It, 
however, turns out that under the assumption of pseudo-Dirac, all probabilities of neutrino 
oscillations are describable in terms of 6 mass eigenvalues and just one 3x3 unitary matrix 



U, "MNS matrix" ||15||, which has three mixing angles and one CP violating phase and just 



corresponds to the CKM matrix in quark sector. 

To see this let us now discuss the diagonalization of the mass matrix M. Ignoring mag- 
netic field, only chirality preserving transitions are important. So what we should diagonalize 
is M^M, rather than M itself. Keeping terms up to the first order in Majorana masses, which 
may be justified in the pseudo-Dirac hypothesis, M^M reads as 

M^M ~ ( M ° Md M* L Mh + MlM* R \ 

~ \ M* D M L + M R M D M* D M l D J' y,) 

Because of pseudo-Dirac property, the dominant matrix is Mjj. Thus we first diagonalize it 
by bi-unitary transformation; 

W R M D U L = diag(rai, m 2 ,m 3 ) = M, 

v olL = {UlJoo. v iLi V a L = {U R ) a i (8) 

Accordingly, M^M is cast into the following form by a unitary transformation due to a 6 x 6 
unitary matrix V, 



. r2 rrt n a\ tt* tC/t i tC/ttt^ 

V j (M f M) V 



M 2 UlMlU* L M + MU R M* R U* R 

MUlM L U L + U R M R U R M ' M 2 

where Ul has been rewritten simply as U, as it is the only matrix which appears in the 
formulae of neutrino oscillations. M 2 is a diagonal matrix, while the matrices in off-diagonal 
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position, e.g. MXJ\M L XJ L + UrMrUrM, have not been diagonalized yet. Such off-diagonal 
matrices, however, are much smaller than M 2 due to the pseudo-Dirac property and seem to 
be negligible, anyway. It is not quite right, since once we ignore these off-diagonal matrices 
there appear degenerate pairs in the eigenvalues of M^M, i.e., each of m 2 (i = 1,2,3) 
appear twice. It is a general wisdom in perturbation theory that when there is a degenaracy 
in eigenvalues, eigenvalues and eigenstates can be fixed only after we include first order 
perturbation, that connects the members of the pair, while other perturbations connecting 
different pairs may be safely ignored. This means that V* (M* M) V can be effectively 
decomposed into three independent block-diagonal matrices. Each block- diagonal matrix 
takes a form of, 



TO- 
TOi€i 



rriiej 

2 

mf 



(2=1,2,3), 



(10) 



where e« = {U^MlUl + U r MrUr)u, and \ei\ <C rro,, because of the pseudo-Dirac property. 
It is now easy to see that we obtain, in total, 6 mass eigenstates 



1 



ViS 



1 



. Their mass-eigenvalues are given as 



ViA 



\y%L 



where e % 



e l ^V iL ), 



1,2,3), 



TO, 



iS 



TO, 



to,- e, ; 



rn 



iA 



to„- 



TO,; ti 



(2 = 1,2,3). 



To summarize, i]) w l is related to mass-eigenstates as, 



wL 



V 



ViS 
ViA 



;i2) 



(13) 
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-i(j> 3 



(14) 



and correspondingly the mass-squared matrix is diagonalized as 

V ] (M f M) V — M 2 ^ M di = diag(TOis, m 2S , m 3S , m 1A , m 2A , m 3A ) 



(15) 



Now the neutrinos emitted by weak interactions (weak eigenstates) are expressed in terms 
of mass eigenstates Uj S , v^ A (j — 1 — 3) and a unitary matrix U, as follows: 



V<xL = U, 



aj 



Vj S + WjA 

V2 ' 



(16) 



The fact that there appears only single 3x3 unitary matrix U, even though we started from 
an arbitrary 6x6 mass matrix M, is one of our main results based on the pseudo-Dirac 
property. 
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3. Formulae for neutrino oscillations 



We will now derive the formulae for neutrino oscillations in terms of the differences of 
6 mass-squared and single unitary matrix U. Though there is no reason to expect apriori 
some specific pattern of neutrino masses, we can still get some useful information on the 
pattern from the reported data on neutrino oscillations fl[], 0, ||. Namely, once we regard 
the mixing angles in U as small, as suggested by the CKM matrix in quark sector, the 
mass-squared difference, responsible for each observed neutrino oscillation, is given as 

solar neutrino : mi|ei| ~ 1CT 5 — 10~ (eV 2 ), 

atmospheric neutrino : m2|e2| ~ 1CT 3 — 10~ 2 (eV 2 ), 

LSND : Am 2 12 ~ KT 1 - l(eV 2 ). (17) 
This knowledge suggests (with a little prejudice) a hierarchical structure of mass differences 

mi|ei| <C m 2 |e 2 | <C m 3 \e-i\ <^ Am 2 2 Am^, (18) 

where Am?- = m 2 — mf. The hierarchical structure makes the formulae for neutrino oscilla- 
tions simple and easy to be compared with the data. 

3.1 A general formula for vacuum oscillation 

We first note that except for the case of solar neutrino oscillation, both of atmospheric 
neutrino oscillation and the oscillation in LSND experiment are well described by vacuum 
oscillations. This is definitely true for LSND case, but may need some care in the case of 
atmospheric neutrino, as will be commented on below. 

In general, the probability of finding a state, born as an active state u a i at time 0, in an 
active state (a and (3 may be the same) at time t is given by 

A//^ 3 2 2 

P(u a - u p ) = \(VeMi^^t}V% a \ 2 = \\Y.U^M^t)+^v{^t)W a3 \ 2 . (19) 

i=i 

3.2 Formulae for atmospheric neutrino oscillation 

As the oscillation of atmospheric neutrino is sensitive to the mass difference m 2 |e 2 | ~ 
10~ 3 (eV 2 ) II, under the mass hierarchy eq.(18) mi|ei| may be ignored and V\ can be 
regarded as pure Dirac particle, i.e. v\ = ^= {{vis + c - c -) + i{viA + c.c.)} with a unique 
mass mis — m iA = mi. The oscillation of atmospheric neutrino is due to the interference 
between z/ 2 s and u^A, an d the matter waves of other states, ui, u-ss and z/ 3j 4, do not interfere 
with U2s, ^2A or with each other, when time-average is taken for the high frequency modes in 
the oscillation probability. Thus the formula for the probability of atmospheric to survive 
till time t, relevant for the zenith angle distribution, simply reads as 

P(y»L - zvJatm = l^il 4 + l^ 2 | 4 cos 2 (^H) + (20) 
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Let us note that there are constant terms l^i] 4 , |£^3 1 4 coming from the time-average of the 
high frequency modes, in sharp contrast to the conventional formula in a simplified 2 states 
system. The following formulae are also relevant for the analysis of atmospheric neutrino 
oscillation: 



i>eL-^ e L) atm = I U el | 4 + | U e2 1 W (^H) + \ \ U e3 \ \ (21) 



P( V - ^L) atm = l^i| 2 |^ii 2 + i^ 2 | 2 |^2|W(^|p*) + ^C/ M 3| 2 |C/e3| 2 . (22) 
3.3 Formula for LSND neutrino oscillation 



12 ~ 



As the neutrino oscillation observed by LSND is sensitive to the mass difference Am 2 
1CT 1 — 1 (eV 2 ) ||, under the hierarchy (18) all neutrino states can be regarded as pure Dirac 
particles, i.e. J7ij|ej| = and m 2 s = m 2 A = m 2 (i = 1,2,3). The LSND neutrino oscillation 
is, therefore, due to the interference between v\ and z/ 2 , and the matter waves of another 
state, 1/3, does not interfere with v\,v 2 , when time-average is taken. Thus the formula for 
the transition probability of LSND neutrino simply reads as 

P(^L^^eL)LSND = IM 2 |^ir+IM 2 |^ 2 | 2 + IM 2 |^ 3 | 2 

f\ 777 ^ j\ 777 ^ 

+ 2Re(U» 1 U: i U*U e2 )cos(—^t) - 21m(U fll U: i U; 2 U e2 )sm(—^t) 



Ami 



H\U, 2 U e2 \ 2 + Re(U; 2 U e2 U, 3 U: 3 )}sin 2 (—^t) 



A 777,^ 

+ 2lm(U; 2 U e2 U, 3 U: 3 )sm(^^t)+2\U, 3 U e3 \ 2 . (23) 

3.4 The matter oscillation of solar neutrino 

The time-evolution of the system in the presence of matter effects are governed by 

dt\u aL J \v aL 1 

where 



= H^Md i&g + V*AV}V^ ^ = diag(a e ,a M ,a T ,0,0,0), (25) 

and the elements of the matrix A denoting the matter effects of u e , u^, v T are given as a e = 
^7=(2N e — N n ), a M = a T = ^|(— iV n ). It is easy to see that for the energy range of solar 
neutrinos and the mass hierarchy of (18), the matter effects inside the sun satisfy 

E a a ~ mi\ei\ m 2 \e 2 \ . . . <C Amj 3 (a = e,/i,r). (26) 

Under the hierarchical structure of mass differences and the matter effects, in the Hamilto- 
nian in the base of the mass eigenstates, + V^AV, all off diagonal matrix elements 
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due to the matter effects can be safely ignored except the ones in 2 x 2 subsystem of v\s and 
i>\a- Thus the heavier states, 2^5, v^a (i = 2, 3) are decoupled from the two states subsystem 
of (z^is, P\a)i an d the time-evolution of the subsystem is governed by a 2 x 2 Hamiltonian 

1 tm \ s o^i ( £ 1^.^(1 J). (27) 



2E\ m{ A ) 2 



=e 



This Hamiltonian can be rewritten in the base of (v\l, Vil) (ignoring a piece proportional to 
a unit matrix) as 

fE; =e |£UV ^sin2#\ 



-sin 26 ^cos 20 



_ (2f 

with Am 2 = 2mi|ei| and 9 = ~. This hamiltonian just corresponds to the one in the MSW 
mechanism J7|, though in our case the mixing angle is maximal and the matter effect has 
been modified into Y7 a=e \U a i\ 2( i a ■ Now the survival probability of solar neutrino (time 
averaged) can be written in a simple form 

P(VeL -> ^L) so l ar = P{PlL ~> ^lL) eff + £ 1^1 4 + ^ We?,\\ (29) 

where the survival probability P{v\l — > vn) in the effective 2 states system of (^il?' 7 il) is 
calculable by use of the Hamiltonian (28). A similar reduction formula based on a hierarchical 
mass structure was obtained in Ref.JIl]] in order to reduce solar neutrino oscillation in 3 
generation model into that in an effective 2 generation model. 

4. Pseudo-Dirac scenario confronted by the data on neutrino oscillations 

4.1 Comparing the formulae with the data 

As the application of the formulae we have derived for the neutrino oscillations, we are 
now going to compare them with the corresponding experimental data. As we have already 
advertised we will see that small generation mixing angles in the unitary matrix U is just 
what we need to explain the pattern of mixing angles observed in neutrino oscillations of our 
interest. 

Thus we first consider the case where generation mixings are small, though the formulae 
we have derived above are applicable for arbitrary generation mixings. Retaining only the 
leading contributions for small generation mixing angles (say di, 8 2 , 63 just as the angles in 
CKM matrix), we get the following formulae, relevant for each neutrino oscillation 

solar neutrino : P{v eL -> ^eL) so \ ar - P{v\l -> v \l)^ (30) 

atmospheric neutrino : 1 - Piv^ -> zvJatm - sin 2 ^^ 2 t), 

P{v eL -> z/ eL ) atm ~ 1, 

P{ynL -> ^ex) a tm - °> ( 31 ) 



LSND : P(^ L ^ z/eL ) LSND ^ 4 |f/ e2 | 2 sin 2 (^|it). (32) 



Am 



AE 
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In the formula for LSND {UezU^l has been neglected compared with {U^U^l, as is suggested 
by the hierarchical mixing angles, 9\ ^> 9 2 #3, in quark sector. In the case of solar neutrino 
oscillation, P(uil — > ^ii)eff * s obtainable from the time evolution governed by the effective 
Hamiltonian Eq.(28), which is very similar to that in the MSW mechanism with maximal 
mixing angle, 7r/4 . To be precise, in the Hamiltonian (28) the matter effect is not y/2 Gp N e 
as in the case of MSW. The difference, however, is not large, as long as the generation mixings 
are small and also because the contribution of neutral current, the term proportional to N n , 
is relatively suppressed compared with that of charged current by a factor ~ 1/12. On the 
other hand, very recently we have heard of the news || that the data of super-Kamiokande 
on solar neutrinos favors the MSW-type solution with large mixing angle. Thus our scenario 
of pseudo-Dirac provides a natural framework to derive the large angle solution suggested by 
the data. In the case of atmospheric neutrino oscillation, the factor in front of sin 2 ( m |^ t) 
is 1 = sin 2 (2 x |), which just corresponds to a vacuum oscillation with maximal mixing, 
strongly suggested by the data on atmospheric neutrinos ffl]. 

In a contrary, in the case of LSND our formula gives that of ordinary 2 generation scheme 
with small generation mixing, if we identify 4|?7 e2 | 2 with sin 2 2#i. This is just consistent with 
the experimental data 0, which says sin 2 2#i < 0.04, when combined with the data from 
BUGEY experiment. Another meaningful constraint on the generation mixing may come 
from the data of CHOOZ experiment |F7|]. As the mass-squared difference which is sensitive 
to the CHOOZ experiment is comparable to that in atmospheric neutrino oscillation, we 
may write down a similar formula to Eq.(20) for the disappearance of T7 e , 

1 - P(u eL -> Z4L)CH00Z = 1 - l^ell 4 - |^e 2 | 4 COS 2 (^|^t) - i|f/ e3 | 4 

^ 2\U e2 \ 2 , (33) 

where the unitarity of U matrix was used and an approximation of \U e ^\ <C \U e 2\ was useful 
to simplify the result. The upper bound on the probability of disappearance from CHOOZ 
puts a bound 

\U e2 \ 2 ~ sin 2 #1 < 0.05, (34) 
which is a little weaker than the upper bound stated above. 

4.2 Problems to be settled 

We have seen that the pseudo-Dirac scenario just provides the favored solutions to the 
solar and the atmospheric neutrino problems with (almost) maximal mixings, suggested 
by the recent data, invoking the oscillations mainly into sterile states, while LSND data is 
naturally explained by ordinary generation mixing between active states with a small mixing. 
We naively expect that neutrino oscillations of solar or atmospheric neutrino into an active 
state and a sterile state cannot be clearly discriminated, since basically these experiments 
are disappearance experiments. Roughly speaking this expectation is certainly true, but 
the data from Super-Kamiokande experiment have reached the precision which is enough 
to distinguish these two cases. The most recent data seem to regard the maximal mixing 
solutions of neutrino oscillations into sterile states with disfavour [|], ||. We would like to 
discuss some possible ways to avoid the difficulty for each case. 
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Atmospheric neutrino 

Though we have neglected the matter effect in the atmospheric neutrino oscillation, 
the matter effect of the Earth becomes non-negligible for higher neutrino energies. It has 
been pointed out that oscillations into v T and a sterile state have different zenith-angle 
dependence, as only in the case of the oscillation into the sterile state the matter effect 
affects the time-evolution of the neutrino states. Comparing with the data, combining the 
analysis of neutral current enriched events, the Super-Kamiokande collaboration claims that 
the oscillation into the sterile state with maximal mixing is regarded with disfavour . The 
possible ways to evade this problem, we can think of, are the following, 
a. When the oscillation to the sterile state is analyzed, simplified 2 states system of (u^, vg) 
(ug denoting a sterile state) is assumed. In the scenario of pseudo-Dirac, however, we have 
6 neutrino states to participate in the oscillation, and the formula for the oscillation, as seen 
in Eq.(20), is different from that in the simple 2 states system, typically having additional 
constant terms (for non- vanishing generation mixings). Including the matter effect of the 
Earth, our formula is modified into 



P{v»L - ^)atm = W,2\ A P(y2L - ^2L) eff + |^i| 4 + ~ IM 4 , (35) 

where the survival probability P{v2L — > ^2i) e ff in the effective 2 states system {y^Li v-il) is 
that for the mass-squared difference 2 m 2 e 2 and the matter effect Y7 a=e \U Q .2\ i Q>a- Thus both 
the depletion rate of atmospheric and the zenith angle dependence should be reanalyzed 
by use of this formula before some definite conclusion is derived. 

b. We just would like to point out that there is a claim that — > vg oscillation with almost 
maximal mixing may not be ruled out, even in the simple 2 states system of (z/ M , vg). It 
has been pointed out that a y 2 analysis of the recent data does not exclude the maximal 
v \i — * v s oscillation solution with any significant confidence level, once various theoretical 
uncertainties and experimental systematic errors are included Jl8 . 

Solar neutrino 

The Super-Kamiokande collaboration claims that the solar neutrino oscillation into the 
sterile state with maximal mixing is not favoured 0. The point is that v e — > u^, v T and u e — > 
vg oscillations give slightly different contributions in the Super-Kamiokande detector, as the 
final active states contribute to the event rate, while the sterile state does not. Accordingly, 
the survival probability of solar neutrino should be relatively higher in the case of oscillation 
into the sterile state. Thus the "Be neutrino problem" in the chlorine experiment becomes 
more severe in the sterile case. The possible ways out of this problem, we can think of, are 
the following. 

a. Again the claim that the solar neutrino oscillation to the sterile state is not favoured is 
based on the analysis assuming a simplified 2 states system of (u e ,ug). In the scenario of 
pseudo-Dirac, however, we have 6 neutrino states to participate in the oscillation, and the 
formula Eq.(29) should be utilized to see whether the oscillation into the sterile state can 
accommodate all data of solar neutrino experiments or not. 

b. In the case of solar neutrino oscillation, the presence of solar magnetic field is potentially 
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important, though we have ignored it in the above discussions. It has been pointed out that 
in the presence of the magnetic field, MSW type oscillation u e i —>■ z7 eL may be followed by a 
RSFP type oscillation V eL — > V eR) ignoring generation mixing ||. The final state, now being 
an active state, contributes to the event rate and may remedy the "Be neutrino problem" . 
c. In the above argument, Majorana masses have been treated as small perturbation, but it 
may not be unnatural to expect that in the first generation the Dirac mass is so small that 
the effect of Majorana masses is relatively enhanced, leading to a relatively small mixing 
between the active and sterile states. 

Cosmological problem 

There is another type of problem, i.e. cosmological problem. The well-known limit of 
number of effective neutrino species during nucleosynthesis puts a stringent bounds on the 



mass-squared differences and mixings of neutrino oscillations into sterile states It, 
however, has been pointed out that once relatively large relic neutrino asymmetry L u , say 
L u > 10~ 4 , is realized such problem can be evaded. For the details of the argument refer to 
PUfl. 



We finally briefly comment on the related issue, i.e. neutrino-less double /9-decay. In our 
pseudo-Dirac scenario, the relevant lepton number violating Majorana mass is given as 

i {U ai f (m iS - ma) ^ - {U ai ) 2 (36) 

which is well below the experimental upper bound under the mass hierarchy (17), (18) 
with €i <C rrii and for small generation mixings. This is basically because the pseudo-Dirac 
neutrinos are almost Dirac particles and the lepton number is only slightly violated by their 
masses. 

5. Theoretical framework for the pseudo-Dirac neutrinos 

Though there have already appeared a few attempts to construct models for the pseudo- 



Dirac neutrinos ||12|| , in this article we instead list up the problems to be resolved before a 
realistic model is constructed, and argue about possible theoretical framework to provide 
natural mechanisms to solve the problems. 

First of all it is worth noticing that the smallness of Majorana masses needed to realize the 



pseudo-Dirac neutrinos satisfies the naturalness condition of 't Hooft ||21[| , since if Majorana 
masses are absent, Ml = Mr = 0, the symmetry of the theory, i.e. the lepton number 
symmetry is enhanced. So the smallness will be stable under the radiative correction. 

We, however, still have the following problems to be settled at the classical level: 

1. How to explain the relation Mr -C Mp ? 

2. How to explain the relation M^ -C Mp ? 

3. How to explain the smallness of Mp itself ? 

We will discuss the possible theoretical frameworks to resolve these problems successively 
below. 

1. Problem 1 
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a. 4-dimensional framework 



In considering the possible theoretical framework to resolve this problem, it may be help- 
ful to reconsider the conventional see-saw mechanism in the language of gauge invariant 
operators. Suppose that the gauge symmetry of our world is U(l) em (QED), then the Majo- 
rana mass term ttil v\ is gauge invariant. Therefore there will be no reason to expect that 
mi should be small. As the matter of fact, the gauge symmetry of the standard model does 
not allow the mass operator, as it is gauge variant. Thus the Majorana mass is provided 
by an irrelevant operator jh L\ H 2 , where H denotes the Higgs doublet and L L = (vlJlY 
is a lepton doublet. As M is a gauge invariant mass, it can be arbitrarily large. Thus the 
essence of see-saw mechanism may be understood as the decoupling of some gauge singlet 
heavy particle with mass M (which need not to be Vr). It, therefore, will be not unnatural 
to expect that a similar thing happens for the right-handed Majorana masses, as well. The 
Majorana mass term uir v\ is gauge invariant in the standard model and rrtR is regarded 
to be quite large. This may not be true in some physics beyond the standard model. For 
instance in the left-right symmetric model SU(2) L x 811(2)^ x U(l) [0] with a SU(2) R 
doublet Higgs H R (in addition to the ordinary SU(2)i doublet Higgs), the m R v R operator 
is no longer gauge invariant, and will be replaced by an irrelevant operator -gj L R H R , with 
Lr = (v R ,l R y. Thus the decoupling of a gauge singlet heavy particle whose mass M' is 
much larger than the scale of SU(2)/j breaking, M' ^> (Hr), may imply the smallness of 
rriR. The heavy particle can be identified with a gauge singlet fermion S, having a Yukawa 
coupling Lr S Hr and a large Majorana mass M'. A diagram with the exchange of S yields 
the irrelevant operator L R H R . 

b. Framework with extra dimensions 

Recently there has appeared revived interest in higher dimensional theories with extra 
dimensions as a possible solution to hierarchy problem [ 23"| , [23 , [25 1 . In the scenario of large 



extra dimension p3| |, the higher dimensional Planck scale is regarded as comparable to the 



weak scale, while in the scenario of small extra dimension [^5j all masses in the visible brane 
are claimed to be strongly suppressed by the "warp factor". Therefore the conventional 
see-saw mechanism ||, which needs large right-handed Majorana masses (tur 3> M^y), may 
not work. In the case of pseudo-Dirac scenario, on the contrary, what we need is very 
small or even vanishing rriR. It is interesting to note that the presence of extra-dimensions 
may provide a natural mechanism to realize this. For instance, let us consider a theory 
in 5-dimensional space-time, where the particles of the standard model , including ul, are 
assumed to reside on a 3-brane, while gauge singlet fields, such as vr, may reside in the 
bulk H26fl , WA - We note that in 5-dimensional space-time Majorana spinor is known not to 
exist. (Majorana spinors exist only in the space-time of D = 2,3,4 (mod 8).) Thus the 
mass term rriR v\ + h.c. is not Lorentz invariant in this space-time and we have vanishing 
tur. Actually in such odd dimensional space-time chiral fermions do not exist, and vr 
should be accompanied by a gauge singlet Dl to form a full spinor ip = (vr, vl) 1 - The full 
spinor may have a large Dirac mass and vr may be decoupled from low energy phenomena. 
This difficulty may be evaded when the extra space is an orbifold S /Z2, the extra-space 



suggested by the recent works [E3|, [55J, PBfl. This is because the discrete symmetry Z 2 
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just corresponds to a symmetry under the transformation ijj — > 75^, which in turn behaves 
as a chiral transformation in the 4-dimensional sense, thus making the Dirac mass term 
prohibited. Strictly speaking, ip may also have masses due to non-zero K-K modes. For 
relatively small sizes of the extra space, such non-zero modes tend to be decoupled from the 
system and do not significantly affect the above argument. 

2. Problem 2 

In some sense, this problem may not be a real challenge; in any viable model the re- 
lation Ml <C Md must be automatically built in, since otherwise custodial symmetry is 
significantly violated by the VEV of a SU(2) triplet representation and Ap = p — 1 gets 
sizable contribution, in contradiction with the data. Furthermore if we have already got 
some mechanism to realize small but non- vanishing rriR in a mechanism discussed above, we 
are satisfied with vanishing mi just as in the Standard Model. In a model with vanishing 
rriR, however, it becomes crucial for the pseudo-Dirac scenario to slightly violate the lepton 
number by small but non- vanishing rriL- It is worth noticing that even in a higher dimen- 
sional model with brane picture the vl is allowed to live only on the brane, and the issue 
concerning mi is essentially 4-dimensional. 

If we wish to get the small at the classical level as the form of a renormalizable 
operator we should introduce a SU(2) L triplet Higgs H T , whose VEV should be small not 
to contradict with the custordial symmetry. Then a marginal operator L 2 l Ht gives small 
rriL. Or we may invoke an irrelevant operator i L\ H 2 , as we discussed above, with a large 
mass scale M 3> (H). The operator may be the result of the exchange of a hypothetical 
gauge singlet S', which has a Yukawa coupling LlS'H and a large Majorana mass M. The 
S' should not be identified with Vr. If we work in the framework of extra dimension with a 
brane, however, there may not be a good reason to assume that the gauge singlet S' resides 
only on the brane. Even if tul is forbidden at the classical level, the small w,l may still 
be produced at the loop level, as far as there is some seed to violate the lepton number. 
The prototype model of this kind may be the model where the presence of a charged SU(2) 
singlet scalar violates the lepton number explicitly |29| . 



3. Problem 3 

In pseudo-Dirac scenario, Dirac masses provide mean masses of neutrinos. Thus crucial 
problem is how to explain the smallness of the neutrino Dirac masses compared with those 
of charged leptons or quarks. 

a. 4-dimensional framework 



In the 4-dimensional framework "Dirac see-saw" mechanism has been put forwarded (12 
in which the marginal operator to give Dirac mass term H Vl vr is forbidden, by a discrete 
symmetry, and the Dirac masses are provided by an irrelevant operator with d > 4, just 
as happens in the conventional see-saw mechanism. Thus if the coefficient is sufficiently 
suppressed by the inverse of some large mass scale, we get small Dirac masses. 

b. Framework with extra dimensions 
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Another intriguing possibility to realize the small neutrino Dirac masses is to invoke the 
presence of extra space f26|, [p7fl . For instance in the 4 + n-dimensional theory with a 3- 
brane, the original higher-dimensional Yukawa coupling f H I7 L u R has a Yukawa coupling 



constant fo, which behaves as 1/JMJ, with Mf being fundamental mass scale of the theory. 

The 4-dimensional Yukawa coupling f± is thus given by ~ 1/ JMJV (V : the volume of 

the extra dimension). The factor 1/yV comes from the overlap of the three fields on the 
brane, and may be suppressed by the largeness of the extra dimension. In fact in the 
ADD model |23j, by use of the relation M 2 pl = Mj +n V, we get a n-independent result 



U ~ (37) 
M v { v ; 

which is ~ 10~ 16 for e.g. Mf ~ ITeV. In this way, a small Dirac mass, ~ f\ Mw ~ Ja Mf, 
is achievable. 

Acknowledgment 

The authors would like to thank W.J. Marciano, Y. Okada, O.L.G. Peres, and O. Yasuda 
for useful and informative discussions. One of the authors (C.S.L.) would like to thank 
the members of the BNL theory group for their hospitality where a part of this work was 
completed. This research was supported in part by the Grant-in-Aid for Scientific Research 
of the Ministry of Education, Science and Culture, No.09246105, No.12047219, No.12640275. 



References 

[1] T. Toshito, Proceedings of the 30th International Conference on High Energy Physics, 
Osaka, Japan, 2000; Y. Fukuda et al., Phys. Rev. Lett. 82 (1999) 2644; Phys. Lett.B467 
(1999) 185. 

[2] Y. Takeuchi, Proceedings of the 30th International Conference on High Energy Physics, 
Osaka, Japan, 2000; Super-Kamiokande Collaboration. Phys. Rev. Lett. 82 (1999) 1810; 
Phys. Rev. Lett. 82 (1999) 2430. 

[3] R.L. Imlay, Proceedings of the 30th International Conference on High Energy Physics, 
Osaka, Japan, 2000; C. Athanassopoulos et al., Phys. Rev. C54 (1996) 2685; Phys. 
Rev.C58 (1998) 2489. 

[4] S. Fukuda et al, Phys. Rev. Lett. 85 (2000) 3999, |hep-ex/ 0009001 . 

[5] L. Wolfenstein, Nucl. Phys., B186(1981) 147. 

[6] M. Kobayashi, C.S. Lim and M.M. Nojiri, Phys. Rev. Lett. 67 (1991) 1685. 

[7] S.P Mikheyev and A.Yu. Smirnov, Sov. J. Nucl. Phys. 42 (1985) 913; L. Wolfenstein, 
Phys. Rev.B17 (1978) 2369. 



14 



C.S. Lim and W.J. Marciano, Phys. Rev. D37 (1988) 1368; E. Kh. Akhmedov, Phys. 
Lett.B213 (1988) 64. 

T. Yanagida, Proceedings of Workshop on Unified Theory and Baryon Number in the 
Universe, ed. by 0. Sawada and A. Sugamoto, KEK, Japan, 1979; M. Gell-Mann, 
P. Ramond and R. Slansky, in Supergravity, ed. by P. van Nieuwenhuizen and D.Z. 
Freedman, North Holland Publ. Co., 1979. 

H. Minakata and H. Nunokawa, Phys. Rev. D45 (1992) 3316; C. Giunti, C.W. Kim 
and U.W. Kim, Phys. Rev. D46 (1992) 3034; For earlier related discussion on the 
pseudo-Dirac neutrinos, see also S.T. Petcov, Phys. Lett. B110 (1982) 245. 

R. Foot and R.R. Volkas, Phys. Rev. Lett. 75 (1995) 4350; A. Geiser, Phys. Lett. B444 
(1998) 358; W. Krolikowski, Acta Phys. Polon. B31 (2000) 663, frep-ph/99 10308. 



D. Chang and O.C.W. Kong, Phys.Lett. B477 (2000) 416, jhep-ph/99l2268| 



C. Giunti, M. C. Gonzalez-Garcia, and C. Pena-Garay, [hep-ph/0001101| ; O. Yasuda, 
[hep-ph/ 0006319 . 



A. Hime and N.A. Jelley, Phys. Lett.B257 (1991) 441; B. Sur et al, Phys. Rev. Lett.66 
(1991) 2444. 

Z. Maki, M. Nakagawa, and S. Sakata, Progr. Theor. Phys.28 (1962) 870. 

C.S. Lim, Resonant Solar Neutrino Oscillation versus Laboratory Neutrino Oscilla- 
tion Experiments, in Proceedings of BNL Neutrino Workshop 1987, BNL report 52079 
(1987). 

Chooz Collaboration, M. Apollonio et al, Phys. Lett. B420 (1998) 397. 
R. Foot, |hep-ph/ 0007065 . 



R. Barbieri and A. Dolgov, Phys. Lett.B237 (1990) 440; K. Enqvist, K. Kainulainen, 
and J. Maalampi, Phys. Lett.B249 (1990) 531. 

R. Foot and R.R. Volkas, Phys. Rev. Lett. 75 (1995) 4350. 

G. 't Hooft, The proceedings of Cargese Summer Inst. (1979) 135. 

J.C. Pati and A. Salam, Phys. Hev.DIO (1974) 275; R.N. Mohapatra and J.C. Pati, 
Phys. Rev.Dll (1975) 566, 2558; G. Senjanovic and R.N. Mohapatra, Phys. Rev.D12 
(1975) 1502. 

N. Arkani-Hamed, S. Dimopoulos, and G. Dvali, Phys. Lett. B429 (1998) 263. 

H. Hatanaka, T. Inami, and C.S. Lim, Mod. Phys. Lett. A13 (1998) 2601. 
L. Randall and R. Sundrum, Phys. Rev. Lett. 83 (1999) 3370. 



15 



[26] N. Arkani-Hamed, S. Dimopoulos, G. Dvali, and J. March- Russell, Neutrino Masses 
from Large Extra Dimensions, |hep-ph/ 98 11448 . 



[27] G. Dvali and A.Yu. Smirnov, Nucl. Phys. B563 (1999) 63, |hep-ph/99042TT 
[28] P. Horava and E. Witten, Nucl. Phys. B460 (1996) 506. 
[29] A. Zee, Phys. Lett.B93 (1980) 389. 



16 



